Surface-enhanced Raman spectroscopy (SERS): nonlocal limitations 
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Giant field enhancement and field singularities are a natural consequence of the commonly employed 
local-response framework. We show that a more general nonlocal treatment of the plasmonic response leads 
to new and possibly fundamental limitations on field enhancement with important consequences for our 
understanding of SERS. The intrinsic length scale of the electron gas serves to smear out assumed field 
singularities, leaving the SERS enhancement factor finite even for geometries with infinitely sharp features. 
For silver nano-groove structures, mimicked by periodic arrays of half-cylinders (up to 120 nm in radius), 
we find no enhancement factors exceeding ten orders of magnitude (10 10 ). © 2013 Optical Society of America 



While the Raman response of (bio-) molecules is in- 
herently weak, nanostructures may be used to tailor 
and tremendously enhance the light-matter interactions. 
This is the key electromagnetic element of surface- 
enhanced Raman spectroscopy (SERS) [1]. In particu- 
lar, metallic nanostructures [2] are known to support 
plasmonic field-enhancement phenomena which are ben- 
eficial for SERS [3]. In many cases, field singularities 
arise in geometries with abrupt changes in the surface 
topography. While such singularities constitute the ba- 
sic electromagnetic mechanism behind SERS, the singu- 
larities are on the other hand an inherent consequence 
of the common local-response approximation (LRA) of 
the plasmons [4]. In this Letter, we relax this approxi- 
mation and allow for nonlocal dynamics of the plasmons. 
To illustrate the consequences we revisit the model ge- 
ometry in 1, initially put forward by Garcia- Vidal and 
Pendry [5] to qualitatively explain the electromagnetic 
origin of the large enhancement factors observed exper- 
imentally. The metallic surface topography is composed 
of a periodic structure of infinitely long metallic half- 





R Vf 



Fig. 1. (a) Groove structure formed by an infinite peri- 
odic array of half-cylindrical nanorods. (b) Cross section 
of the unit cell, (c) and (d) Typical electric-field intensity 
and charge distributions for a dipole mode. 



cylinders of radius R, resting shoulder-by-shoulder on a 
semi-infinite metal film. The steep trenches or grooves 
support localized-surface plasmon resonances (LSPR). 
Near the bottom of the groove the surfaces of the two 
touching half-cylinders become tangential to each other 
and a field singularity forms within the traditional LRA 
of the dielectric function. In the common treatment, the 
field enhancement thus eventually turns infinite [6] while 
it remains finite, albeit large, in any experiment reported 
so far. Geometrical smoothening is known to remove 
the singularity within the LRA and in quantitative nu- 
merical studies a rounding needs to be added to make 
numerical convergence feasible [7,8]. Thus, within the 
LRA framework the field enhancement would just grow 
without bound the sharper one could make the geome- 
try confining the plasmon oscillations. Nonlocal effects 
have been shown to result in large blueshifts and con- 
siderably reduced field enhancements (as compared to 
a local description) in conical tips [9], metallic dimers 
involving small gaps below a few nanometers [10, 11], 
or even vanishing gaps [12]. What is the limit in field 
enhancements that can be achieved with (geometrically) 
ideal structures? This question is important not only 
from the fundamental but also from applied perspective, 
as the answer to it would allow one to determine tech- 
nological tolerances in fabrication of nanostructures de- 
signed for achieving record-high field enhancements. In 
this Letter we show how nonlocal response introduces 
a new intrinsic length scale that serves to remove the 
field singularities, leaving field enhancements finite even 
in geometries with arbitrarily sharp changes in the sur- 
face topography. For the particular geometry of Fig. 1 we 
evaluate 7(r,cj) = |E(r,cj)| / |Eo(o;)| and find no (sur- 
face averaged) SERS enhancement factors (7) exceeding 
ten orders of magnitude. 

The electromagnetic response of a metal is commonly 
divided into intraband contributions [13] and the disper- 
sive Drude free-electron response £d(u) — 1 + i-f^ = 
. .2 

, where a is the complex conductivity also 
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Fig. 2. Surface- averaged SERS enhancement factor (7) 
for the case of R = 75 nm with r = 0.1 nm (upper curves) 
and r = 5 nm (lower curves). For comparison, the dashed 
lines show the results of the commonly employed local- 
response approximation. 



appearing in Ohm's law J = crE. We relax the latter 
local-response constitutive equation and turn to a lin- 
earized hydrodynamic nonlocal treatment [10, 11, 14, 15] 
where the usual Maxwell wave equation is coupled to 
a hydrodynamic equation for the current density, see 
Ref. [11] for the full details of our numerical approach. 
This is the simplest non-trivial extension of the common 
LRA Drude model, which in addition to the usual metal 
parameters (lj p , td, etc.) now also carries information 
about the kinetics of the charge carriers at the Fermi 
level. The strength of the nonlocal correction to Ohm's 
law depends on the Fermi velocity vp which introduces 
a new length scale, being a factor vp/c of the free-space 
wavelength A = 2ttc/uj. For the noble metals, vp/c is of 
the order 10 -2 which explains the overall success of the 
LRA. However, when exploiting plasmonics at the true 
nanoscale, effects due to the nonlocal dynamics start to 
manifest themselves. Field-enhancement structures turn 
out to be prime examples of this. 

We consider the metallic groove structure shown in 
Fig. 1 which has previously been considered as a model 
system to mimic corrugated metal surfaces [5]. Alterna- 
tively, it may be viewed as a model for arrays of the 
more recent groove or channel waveguides [7, 16]. In our 
numerical study, the structure is excited by an incom- 
ing plane wave Eo(cj), normal to the substrate and with 
the field polarized perpendicularly to the axis of the 
half-cylinders, i.e. across the groove cross section. No- 
ble metals are common choices for plasmonics and in 
the following we focus our attention on silver [13]. The 
grooves have been shown to support LSPRs [7] which we 
have previously explored in the context of SERS, using 
a LRA and with the necessary addition of geometrical 
smoothening [8]. To quantify the SERS effect and the 
consequences of nanoscale spatial dispersion, we solve 
the nonlocal wave equation numerically [11]. As an ex- 
ample of our results, Fig. 2 shows the spectral depen- 



dence of (7) throughout the visible regime for groove 
structures with R = 75 nm and with a radius of cur- 
vature of the crevice given by r = 0.1 nm. The LSPR 
at A = 700 nm allows the (surface-averaged) Raman 
rate to be enhanced by a factor of 10 8 . For comparison, 
the dashed line shows results when treating the plas- 
monic response within the common LRA. In both cases, 
the resonant behavior is well pronounced, being caused 
by interference of the incoming field with the gap sur- 
face plasmon mode reflected at the bottom, similarly to 
that described for V-grooves [16]. As a general finger- 
print of nonlocal response, the peak is blueshifted com- 
pared to the expectations from a local-response treat- 
ment of the problem (this happens due to a decrease in 
the gap plasmon index caused by nonlocal effects [10]). 
In this particular case, the LSPR by the common treat- 
ment is off by more than 25 nm which illustrates the im- 
portance of nonlocal effects for quantitative SERS pre- 
dictions. Even more importantly, the common LRA is 
seen to significantly overestimate the enhancement fac- 
tor; for some wavelengths by more than one order of 
magnitude. The large quantitative differences between 
the nonlocal treatment and the traditional LRA are as- 
sociated with changes in the induced-charge distribution 
(insets of Fig. 2). In the common treatment, the charge 
is strictly a surface charge while in the general nonlocal 
case the intrinsic scale vp/uj serves to spatially smear 
out the charge distribution. Effectively, this smearing 
increases the electric field penetration into metal (sil- 
ver) and thereby increases the field absorption (ohmic 
loss) and damping of resonant oscillations. Interpreting 
the field enhancement in a capacitor picture, the finite 
thickness of the charge distribution near the surface in- 
creases the effective separation (beyond that given by 
the metal-surface geometry) and consequently the ca- 
pacitor supports a lower electrical field compared to in 
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Fig. 3. Surface- aver aged SERS enhancement factor (7) 
for the case of R = 15 nm and with r varying from 
1 nm to nm. The dashed line connecting fundamental 
dipole resonances for different values of r serves as a 
guide to the eyes, clearly illustrating both a redshift and 
the saturation effect in the field enhancement as r — » 0. 
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Fig. 4. Near-resonance plots of the surface- averaged 
SERS enhancement factor (7) for arbitrarily well-defined 
grooves without smoothening (r = 0) for six cases with 
R varying from 30 to 120 nm. The inset shows the field- 
amplitude distribution |E|/|Eq| for R = 75 nm. 



the LRA. In general, the intrinsic length scale of the elec- 
tron gas allows to resolve the field also in the proximity 
of very sharp corners and tips. On the other hand, by 
relaxing the sharpness of the trench the influence of spa- 
tial dispersion becomes less pronounced, as illustrated 
in Fig. 2 in the lower set of curves (r = 5nm) where 
the LRA accounts well for the results obtained from a 
full nonlocal treatment. We note a drastic change in the 
field enhancement spectrum, with the fundamental res- 
onance now appearing at around 450 nm, due to a very 
rapid decrease in the gap plasmon index when the gap 
width increases (at the groove bottom) from 0.1 to 5 nm. 

With less geometrical smoothening (i.e. when r is 
made smaller and smaller) the shortcomings of the 
LRA become more severe. The LRA anticipates a 
monotonously increasing enhancement factor [8] and 
decreasing r also causes a stronger interaction be- 
tween neighboring half-cylinders and consequently a red- 
shift [5]. Note that in the interpretation based on gap 
surface plasmons [16], the redshift is simply related to an 
increase in the gap plasmon index when the gap width 
decreases at the groove bottom. In Fig. 3 we decrease 
r from 1 nm down to zero and see how nonlocal effects 
cause a different trend (indicated by the dashed line) 
due to the competing length scales. In particular, for 
r ;$ vf/oj there is a fundamental saturation of the en- 
hancement factor rather than a monotonous increase and 
for our particular choice of the cylinder radius R we see 
that the (7) does not exceed 2 x 10 9 . 

To explore the ultimate limitations on the SERS in 
this geometry, Fig. 4 shows results where we have com- 
pletely refrained from any geometrical smoothening (r = 
0) and where vp/oo is the only length scale that puts fun- 
damental limitations on the field enhancement. As the 
radius R of the half-cylinders is increased from 30 nm to 
120 nm we see a redshift of the peak as also anticipated 
in the LRA [8]. At the same time, the enhancement fac- 



tor exhibits an increasing trend where larger cylinders 
support larger field enhancement by harvesting the in- 
coming field from larger areas. We emphasize that in all 
examples the field enhancement remains finite despite 
the fact that the crevice is arbitrarily sharp and well de- 
fined (r = 0). For the largest radius R considered the 
electromagnetic SERS enhancement factor does not ex- 
ceed 2 x 10 10 . This illustrates the fundamental limita- 
tions imposed by nonlocal response in our specific SERS 
configuration. 

In conclusion, we have shown that a nonlocal treat- 
ment of the plasmonic response leads to new and possibly 
fundamental limitations on the electromagnetic SERS 
enhancement factor, thereby completely changing the 
message of the commonly employed local-response ap- 
proximation of the plasmons. The intrinsic length scale of 
the electron gas serves to smear out the field singularity 
that otherwise would arise from a local-response treat- 
ment and as a consequence the enhancement remains 
finite even for geometries with infinitely sharp features. 
Finally, beyond the linear response fundamental limita- 
tions may arise due to nonlinearities [17]. 
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